Secondary Neutrinos from Tau Neutrino Interactions in Earth 
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The energy dependence of "secondary" neutrinos from the process (v T — > r 
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tau neutrino fluxes (F„ ~ E^ 1 and E„ 2 ), assumed to have been produced via neutrino oscillations 
from extragalactic sources, is evaluated to assess the impact of secondary neutrinos on upward 
muon rates in a km 3 detector. We show that the secondary fluxes are considerably suppressed 
for the steeper flux, and even for fluxes ~ E^ 1 , the secondary flux will be difficult to observe 
experimentally. 



Evidence of neutrino oscillations from measurements of 
the atmospheric v e and fluxes leads one to the conclu- 
sion that v ^ — > v T over distances characterized by the ra- 
dius of the Earth for neutrino energies ~ 1 GeV §. The 
results of a two-flavor analysis yield Am 2 = 2.5 x 10~ 3 
eV 2 and bi-maximal mixing, sin 2 26 = 1. Three flavor 
analyses are consistent with this result ||. 

Atmospheric neutrinos come from cosmic ray interac- 
tions with air nuclei, yielding hadrons, especially mesons, 
which decay to neutrinos. Neutrinos may also be pro- 
duced at the sources of cosmic rays, where energetic pro- 
tons interact with nucleons and photons at the source. 
The pions and kaons produced at these distant astro- 
physical sources are parents of neutrino fluxes. Active 
galactic nuclei and gamma ray bursters are two proposed 
astrophysical neutrino sources j| . Other possible sources 
include exotic particle annihilations Q . Bi-maximal mix- 
ing, in the context of astrophysical neutrino sources, re- 
sults in flavor oscillation from the source ratio of fluxes: 
v e : : v T = 1 : 2 : to flux ratios at the Earth of 
i/ e : : v T = X : 1 : X, independent of the neutrino 
energy, given that astronomical distances are so large || . 

Halzen and Saltzberg pointed out in Ref. Q that high 
energy v T flux attenuation in the Earth differs from 
flux attenuation due to the fact that the r produced in 
charged-current (CC) interactions with nucleons decays 
before it loses energy. For each v T lost in CC interactions, 
another v T appears following each r decay, albeit at a 
lower energy. Detailed evaluations of v T flux attenuation 
in the Earth appear in Refs. [^-^). Depending on the 
incident flux, the v T flux shows a degree of "pile-up" as 
neutrinos of sufficiently high energy interact in the Earth 
and yield neutrinos at lower energies. 

In a recent paper, Beacom, Crotty and Kolb Jl(J] have 
suggested that in addition to a pile-up of tau neutrinos, 
the signal of astrophysical tau neutrinos will be enhanced 
by the appearance of "secondary" neutrinos. These sec- 
ondary neutrinos come from purely leptonic decays of r's. 
The idea is that while i/g {I = e, fi) fluxes starting, for 
example, at nadir angle are extinguished for sufficiently 
high energies, they will be regenerated by the B = 0.18 
branching fraction for r — ► v T £i't, the r being produced 
by v T — > r CC interactions. The flux of vt is less likely to 



be extinguished due to the shorter path-length through 
the Earth (the v T already had to travel its "interaction 
distance,") and its lower energy due to the combined en- 
ergy loss in the CC process and the decay of the r. 

The evaluations in Ref. [fiof were for mono-energetic 
neutrinos. In this paper, we consider a few energy de- 
pendences for the v r fluxes. We show that the sec- 
ondary fluxes are considerably suppressed for steeply 
falling fluxes and even for fluxes ~ E„, the secondary 
flux will be difficult to observe experimentally. 

Neutrino attenuation in the Earth is governed by a 
coupled set of partial differential equations. To illustrate, 
we write the coupled equations for the v T flux. For energy 
dependent flux F Vr (neutrinos/(cm 2 s sr GeV)), 
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Here, £™* = I/Na&uN JllJ an d similarly for the r inter- 
action length, and £^ ec = ^crp for density p and Lorentz 
factor 7 = E T /m T . The quantity X is the column depth 
(in g/cm 2 ), and for example, 
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The cross section normalized energy distribution of neu- 
trinos with incident energy E y and final energy E is rep- 
resented by dn NC /dE in Eq. (2). 

There are similar equations for v T and f fluxes. In our 
discussion of flux ratios in the introduction, we didn't 
distinguish between neutrinos and anti-neutrinos. Mod- 
els predict that F v — F^ to a good approximation. For 
our discussion in this section, we distinguish between par- 
ticles and antiparticles since the secondary neutrinos are 
actually anti- neutrinos (pg) for incident v T . One must 
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keep in mind that there are equal incident neutrino and 
antineutrino fluxes of all three flavors in the context of 
neutrino oscillations. In the event rates evaluated below, 
we sum both neutrino and antineutrino contributions. 

For nadir angle 0, C 1 ™ 1 equals the Earth's diameter 
when E u ~ 40 TeV. Tau neutrino pile-up and associ- 
ated secondary anti-neutrino production is relevant at 
this energy and higher. At very high energies, we have 
shown that the effect of attenuation for a surface detec- 
tor viewing at zero nadir angle is significant even for v T . 
For example, for a flux (with a smooth cutoff at 

E u = 10 s GeV), the attenuated flux of v T at E v = 10 6 
GeV is only 6% of the incident flux at that same energy 
0. The u e and fluxes are about 1% of the incident 
flux at the same energy. Steeper fluxes have an even more 
significant attenuation and thus will have lower relative 
detection rates at high energies, whether it be v e , or 
v T incident fluxes. 

Guided by the falling fluxes and the increased atten- 
uation, wc confine our attention to the energy range 
of E v = 10 3 - 10 8 GeV. For energies below 10 s GeV, 
£int > C dec starting at E ~ 10 8 GeV, G r * T becomes 
important jL2|, however, we neglect it here and confine 
our attention to lower energies. We also neglect the con- 
tribution of C" lt in Eq. (1). With the fluxes considered 
below, the error in this approximation should be small 
for rates evaluated with minimum energies of 10 4 — 10° 
GeV. With these approximations, the set of differential 
equations is simplified and solved [Q using a modification 
of the iterative method detailed by Naumov and Perrone 
in Ref. |l3). For fluxes, G M ^ M from electromagnetic 
muon energy loss effectively eliminates any return of v ^ 
from CC interactions followed by [i decay in the energy 
range of interest. 

The t flux solution to Eq. (1) is responsible for gen- 
crating the secondary neutrino flux. In the energy range 
of interest, one can write the differential equation for the 
i>£ flux as: 

dF Pt (E,X) = F, t (E,X) 

8X £>J U {E) 1 ' 

+ / dE y [G p ^ Pl (E,E v ,X) + G T ^ p '(E,E v ,X)} . 
Je 

In what follows, we set G v ^ v = G p ^ p . At sufficiently 
high energies, neutrino and antineutrino interaction rates 
are equal because the cross sections are dominated by the 
sea quark distributions. The energy at which the interac- 
tion length equals the column depth increases with nadir 
angle, so the approximation is best at larger nadir angles. 
Because the pile-up comes from higher energy neutrino 
interactions, this approximation is not unreasonable, as 
discussed below. 

A second approximation to obtain the secondary flux 
is to take: 

QT^t ^ Q . QT^Vr j ( 4 ) 



for B = 0.18, the branching fraction for r — > I. In fact, 
the \>i spectrum from r decay is a little softer than the v r 
spectrum, so this is an approximation that will slightly 
overestimate the secondary flux. 

With these two approximations, one finds that the 
combination A(E, X) = F Vt (E, X) - Fp t {E, X)/B satis- 
fies the same transport equation as F v (E, X). At X = 0, 
A(E, 0) = F Vt (E, 0) = F Vft (E, 0), so we can write 

F- Vt (E, X) ~ B ■ [F Vr (E, X) - F„„ (E, X)) . (5) 

Eq. (5) will be used in what follows to approximate 
the secondary antineutrino flux for two different incident 
spectra: E^ 1 and E~ 2 . 

In Fig. 1, we show the ratio of the attenuated flux to 
the incident flux at nadir angle 9 — 0, 30 and 60 degrees 
for 

F° = F y (E, X ~ 0) ~ Ni/E • 1/(1 + E/Eo) 2 (6) 

where Eq = 10 8 GeV and Ai is a normalization fac- 
tor, and for F® oc 1/E 2 . The dashed curve is the v T 
flux, the dotted curve shows the attenuation of the other 
two neutrino species. The solid curve, the result of Eq. 
(5), is the secondary T>g flux, were £ — e or /i. Except 
for the smallest nadir angles for the E~ x spectrum, the 
secondary antineutrino flux is a small correction to the 
primary attenuated electron neutrino or muon neutrino 
flux. Even at nadir angle zero, the secondary flux is neg- 
ligible compared to the transmitted primary flux for the 
1/E 2 . spectrum. 

One should note that even though the v T flux domi- 
nates the primary and secondary fluxes, it does not 
dominate the contributions to the muon event rate be- 
cause of the branching fraction B — 0.18 of r — ► n to- 
gether with the effect of energy loss as the v T converts to 
a t which then decays to a p. 

As a quantitative illustration of the implications of 
the secondary neutrino flux, we consider the muon event 
rate from the F° +[ , ~ N\/E v case (Eq. (6)) with 
N\ = 10~ 13 /(GeV cm 2 s sr) for each neutrino flavor. 
The normalization factor Ai is chosen to be in line with 
the Waxman and Bahcall gamma ray burster flux of Ref. 
fl4j| , in which Ai = 4 x 10 -13 (in the same units) for the 
sum of all neutrino species and for E v < 10 5 GeV. The 
event rates shown below are for an underground detector 
of 1 km 2 effective area, for example, the proposed Ice- 
Cube detector [^5). Details of the calculation appear in 
Ref. Our integrals over neutrino energies were per- 
formed from the minimum neutrino energy E^ > 10 4 or 
10 5 GeV, up to a maximum neutrino energy of 10 8 GeV. 

In Fig. 2, we include the contributions from — > fi 
and v T — ► t — ► [.i (shown with the dashed lines) and 
in addition, the corresponding antineutrino induced an- 
timuons from v T — * r — > — > p, (solid line). Anti- 
neutrino induced muons and antimuons are also included. 
We note that for the — > fi rates, doing the neutrinos 
and antineutrinos separately results in at most a ~ 10% 
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FIG. 1. The ratio of the attenuated neutrino flux for v T (dashed), — v £ (dotted) and secondary vt (solid) to incident flux 
F° oc E- 1 (Eq. (6)) and F? tx E~ 2 . 
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correction to the event rates at the energies shown here. 
This small correction is an indication that our approxi- 
mations to obtain the secondary flux of neutrinos are not 
unreasonable. 
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2. The muon event rate for muons with energy above 
— 10 4 and 10 5 GeV originating from neutrinos with 
E v < 10 8 GeV from M -» fi) + (y T -> r -» fi) plus 
antiparticles (dashed) and additionally from secondaries via 
(y T — > t — > — y /i) plus antiparticles (solid) for the incident 
flux as in Eq. (6) and Ni = 10 _13 /(GeV cm 2 s sr) for each 
incident neutrino plus antineutrino flavor. 

The secondary neutrino contribution to the muon 
event rate has its largest relative contribution at nadir 
angle zero, with an enhancement over the — > [i plus 
► t ► fi rate of 50-60% for the \/E v flux. At this 
angle, the ratio of rates of the secondary contribution to 
the muon event rate relative to the event rate of muons 
from tau decays is quite large, about 1.5 for > 10 4 
GeV and 2.6 for > 10 5 GeV. Unfortunately, this is 
where the event rate is smallest and statistics are low. By 
a nadir angle of ~ 60° (1 rad), where the event rate is 
roughly a factor of 10-20 larger, depending on the mini- 
mum muon energy, the enhancement in the overall muon 
rate is about 25%. At this nadir angle, the secondary v 
produced muons are equal to the tau decay muon rate for 
E, 
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For more rapidly falling fluxes, the contribution of sec- 
ondary neutrinos to the event rates is smaller. The atten- 
uated fluxes shown in Fig. 1 for l/E^ yield secondary en- 
hancements which are quite small. The secondary muon 
rate is only about 10-15% of the primary — > rate 
at nadir angle zero. Typical theoretical neutrino fluxes 
have spectra that lie somewhere between the \/E v and 
1 1 El cases in this energy range [|| . 

In summary, the energy dependence of the incident tau 
neutrino flux is crucial in evaluations of the implications 
of v T interactions to regenerate v T and secondary and 
v e . Neutrino fluxes are attenuated due to their passage 
through the Earth, even the tau neutrinos, thus moder- 
ating tau neutrino contributions to secondary neutrino 
fluxes. Our evaluation has relied on approximations in- 
cluding Eq. (4), setting = and G v ^ v = G^ . If 
muon rates are determined to be large, even near nadir 
angle zero, then a more detailed evaluation of the sec- 
ondary flux may be in order. With our current theo- 
retical expectations for flux normalizations, however, the 
secondary neutrinos coming from r decays will be dim- 
cult to observe experimentally, as they contribute signifi- 
cantly to a muon excess only at small nadir angles where 
the fluxes are already strongly attenuated and for spectra 
like l/E v . 
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